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The extraction algorithms for pulse amplitude and smoothing of energy spectrum have a great influence on 
energy spectrum of y-rays during the digital detection and analysis procedure. For a CdZnTe digital y detector 
system, different extraction algorithms for pulse amplitude and smoothing of energy spectrum are discussed 
in this paper. The results show that extraction of pulse amplitude using the first-order derivative method and 
smoothing of energy spectrum using the wavelet transformation method may obtain energy spectrum with good 


performance. 
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I. INTRODUCTION 


Cadmium zinc telluride (CdZnTe) digital y detector system 
is composed of a CdZnTe detector, a high speed data acqui- 
sition card and relevant software [1, 2]. The data acquisition 
card receives analog signals from the detector and converts 
them into digital signals [3]. The software extracts the sig- 
nal amplitude, analyzes the statistics, and outputs the energy 
spectrum. The energy spectrum depends on the extraction 
algorithm and statistical characteristics of the physical proce- 
dures [4]. For obtaining desired CdZnTe y-spectra, different 
algorithms to extract the pulse amplitude and to smoothing 
the energy spectrum, and the algorithm optimization, are in- 
vestigated in this paper. 


Il. METHODOLOGY 


The method to obtain the pulse amplitude affects the peak 
position and energy resolution of the spectrum. We studied 
three amplitude extraction algorithms to find a suitable one. 


A. Amplitude extraction algorithms 
1. Extremum method (Fig. 1(a)) 


In this method, each pulse maximums ymax and pulse mini- 
mum Ymin are extracted, and A = 0.9(Ymax — Ymin) is the pulse 
amplitude, considering the impact of noises and backflushing. 
The energy spectrum is obtained by statistically treating all 
the pulse amplitudes. 


2. First derivative method (Fig. I(b)) 


Using the first derivative of pulse amplitudes, channel ad- 
dresses of the maximum values can be achieved [5]. The 
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channel addresses of tıs and tip in Fig. 1(a) corresponding to 
zero value of the first derivative are the initial and final chan- 
nel, respectively. The amplitude values A can be extracted 
from the addresses of (tis, Yıs) and (tib, Yip) Of the original 
pulses: A = yis — Yb. 


3. Second derivative method (Fig. I(c)) 


The channel addresses of the maximum and minimum val- 
ues of the second derivative amplitudes are obtained. The 
initial and final channel addresses of t2, and t» in Fig. 1(a) 
are acquired. The amplitude values A can be extracted from 
the addresses of (tas, y25) and (tap, Yon) of the original pulses: 
A = Ys — Y». 
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Fig. 1. (Color online) Illustration of the methods. (a) original pulse, 
(b) the first derivative pulse, (c) the second derivative pulse. 


B. Filtering algorithm for energy spectrum 
1. Wavelet transformation 


Wavelet transformation is applied to decompose and re- 
construct the pulse signals [6]. To suppress signal noises in 
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the reconstruction process, coefficients of the branches with 
larger noise are reduced, while coefficients of the branches 
with larger signal feature are kept [7]. The processes of sup- 
pressing noise are as follows [8—10]: 


(1) The sym4 wavelet basis function is applied to decompose 
the original signals into three parts. 


rO) = $` ho(n — 2k)ait1 
dP = S h(n — 2k)ai*! 
h(n) = (—1)"ho(N — n) 


(§20,j¢Z), W 


where xË ) and dË ) are respectively the coefficients of 
low and high frequencies under the decomposition scale 
of j. ho(n) and hi(n) are the coefficients of low-pass 
filter and high-pass filter. 


(2) Soft threshold method is applied to contract the high fre- 
quency coefficients of the three parts. 


(3) According to the low frequency coefficients for the third 
parts and the high frequency coefficients for the three 
parts treated by the threshold values, the wavelet recon- 
struction is performed to restore the original signal. 


Wavelet reconstruction formula is given in Eq. (2): 


at! — Y ho(n — 2k)a? +S hlin kad. (2) 
k k 


2. Finite impulse response (FIR) filter 


The basic idea of FIR filter is separating the low frequency 
signals from high frequency noises with an appropriate win- 
dow function and cut-off frequency [11]. With an ideal 
filter frequency response H,4(e)”), an FIR filter frequency re- 
sponse Hy(e?”) = h(n)e~J* (n = 0 > N — 1) is applied to 
approach Hy(e%”). A limited sequence of window function 
w(n) is used to cut off ha(n): 


h(n) = w(n)ha(n). (3) 


In this paper, the sampling frequency is f =10° Hz, with 
the angular frequency of Q, = 27 x 108 rad/sec. The cut- 
off frequency of transmission band is f =10° Hz, with the 
angular frequency being Qp = 27 x 10% rad/sec. The initial 
frequency of the stop band is f = 210° Hz, with the angular 
frequency of Qs = 27x 10? rad/sec. The attenuation constant 
of the stop band is not less than —50 dB. 

For the hamming window in the data processing, the min- 
imum value of decay constant is —53 dB, and the stop-band 
is larger than —50 dB, so the hamming window function is 
chosen as 


w(n) = fossa — 0.46 Ea] Ry(n). A 
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The amplitude function of frequency response is: 
W (n) = 0.54Wpr(w) 


27 2T 
+0.23 [wa (w- x3) We (0+55)| 
~ 0.54Wp(w) 


27 2T 
2 — — — JI. 
+0 3 [wa (2 7) +We (w+ FE) 


Using this method, 99.6% energy can be focused in the 
main lobe of the window spectrum. The width of main lobe 
is 27t/N. And the peak value of the side lobe is less than 1% 
of the peak value of main lobe. 


3. Least square smooth 


The least square smooth method is as follows. To obtain the 
m point in the smoothed spectrum, k + 1 points are selected 
on each side of the m'" point to form a window including 
2k + 1 data points. In this window, the polynomial fitting 
method is applied to fit the original data, and the value of m" 
point of the fitting polynomial is the value of the smoothed 
spectrum at the m" point. As the m point moves, the entire 
smoothed spectrum is obtained [12]. 

Let ym be the original data and Ym be the smoothed data, 
the smoothed value of the m® point is 


Ym = S(2)|e=m = 40, (6) 

as the q order polynomial is used to approach ym 
S(x) = agt+ai(a—m)+a9(a4—m)?+-+-+a4g(x—m)!. (7) 
We use 11 points smoothing (k = 5, a window of 11 data 
points). Quadratic polynomial is applied to fit the data, and 


the least square method is used to obtain the coefficients of 
the polynomial. The 11 points smoothing formula is 


Ym = ( = 36Ym—5 Se 9Ym—A + 44Ym-3 + 69Yym—2 
+ 84Ym—1 + 89Ym + 84Ym+1 + 69Ym+2 (8) 
+ 44ym+3 + 9Ym44 — 36Ym45)/429. 


HI. COMPARISON OF ALGORITHMS 


CdZnTe digital y detector system is used to measure °°Co 
y-rays, and algorithms are compared and assessed for their 
capabilities of obtaining the energy spectra of different am- 
plitudes and smoothing results [13]. The above algorithms 
are applied to treat the obtained impulse data of y-rays. Sev- 
eral metrics are used to evaluate the effects of different algo- 
rithms such as energy resolution, integral nonlinearity peak- 
to-Compton ratio, peak area, etc. 
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TABLE 1. Parameters calculated with the 1“ derivative (FDM), the 2"4 derivative (SDM) and the extreme methods (EM) 


Methods Energy Peak height Peak channel FWHM Energy resolution Peak area Peak-to-Compton 
(MeV) (%) ratio 
FDM 1.17 71.1 319 11 3.44 177.7 1.7 
1.33 49.7 364 13 3.56 97.6 
SDM 1.17 60.7 320 15 4.67 124.5 1.6 
1.33 46.6 366 15 4.10 120.8 
EM 1.17 49.3 318 27 8.52 54.4 1.1 
1.33 32.2 362 20 5.54 29.6 


*All units are dimensionless. 


TABLE 2. Smooth results by wavelet transform (WT), filter smooth (FS), least square moving smooth (LSMS), for 6°Co y-rays 


Methods e, Peak height Epal k a talo Peak channel 

Original 1.17 82 319 
1.33 57 364 

WT 1.17 71.1 10.8 319 
1.33 49.7 73 364 

FS 1.17 69.0 11.0 322 
1.33 47.8 9.2 368 

LSMS 1.17 69.0 11.0 318 
1.33 48.3 77 365 
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Fig. 2. (Color online) Energy spectra obtained with (a) the first 
derivative method, (b) the second derivative and (c) the extremum 
methods. 


A. Acquirement and evaluation of impulse amplitude 


The energy spectrum obtained is relatively rough. In or- 
der to obtain the typical parameters and reduce the influence 
of statistical error on computation of the parameters, wavelet 
smoothing is performed on the spectrum [14, 15]. The spec- 
tra obtained with the first derivative, the second derivative and 
the extremum methods are shown in Fig. 2. 

The spectrum obtained by the first derivative method 
(Fig. 2(a)) not only is smooth, but also has prominent peaks. 
The spectrum obtained by the second derivative method 
(Fig. 2(b)) is distorted on the Compton plateau region and on 
the peak position. The spectrum obtained by the extremum 
method (Fig. 2(c)) is similar to Fig. 2(a), but the peak value is 


Channel drift 
ratio (%) ratio 


Energy resolution Peak-to-Compton 


FWHM Peak area 


0 11 3.44 177.7 1.7 
0 13 3.56 97.6 
3 11 3.41 151.5 1.5 
4 12 3.26 106.0 
1 12 3.76 166.1 1.6 
1 13 3:95 120.1 
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Fig. 3. (Color online) Original energy spectrum(a), and energy spec- 
tra smoothed with (b) the wavelet, (c) the filter and (d) the least 
square methods. 


cut down. The parameters calculated with the three methods 
are listed in Table 1. 

From Table 1, the spectrum obtained by the first derivative 
method is the best of all methods, with the best energy resolu- 
tion, the highest peak height, the maximal peak-to-Compton 
ratio and the maximal peak area. 


B. Smoothness and evaluation of energy spectrum 


The energy spectrum smoothed with the three methods and 
the unprocessed spectra are shown in Fig. 3. 
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The spectrum smoothed by the wavelet method (Fig. 3(b)) 
not only is smooth, but also keeps the feature of the origi- 
nal spectrum without peak position shift. The spectrum pro- 
cessed by the filter method (Fig. 3(c)) is smooth, with shifted 
peak position, though. The smoothed spectrum with the least 
square method (Fig. 3(c)) is not good enough because of some 
glitch impulses. The parameters of the smoothed spectrum 
are listed in Table 2. 


From Table 2, statistical fluctuation of the spectrum can- 
not be removed, or original shape of the spectrum is not kept, 
using the least square smooth method, due to limitation of 
the base vector of the discrete Fourier transform. The inte- 
gration nonlinearity cannot be suppressed effectively by FIR 
filter, the peak position differs remarkably from the original 
position. Because of the multi-resolution for the base vector 
of wavelet transform method, high resolution can be achieved 
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in channel domain and frequency domain, and features of the 
slow and rapid changes also can be fully extracted, hence an 
effective method for smoothing y-ray spectra. 


IV. CONCLUSION 


In this paper, integration nonlinearity, peak area, peak-to- 
Compton ratio are set as the standards of evaluating the exper- 
imental data on the extraction algorithms for impulse ampli- 
tude and the smoothing algorithms for the energy spectrum 
obtained by CdZnTe digital y-ray detector system. The re- 
sults show that using the first derivative method to extract 
the impulse amplitude and the wavelet transform method to 
smooth the energy spectrum have obtained the better results 
taking no account of the computational time. 
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